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Abstract We investigate the Mbh — c* relation for radio-loud quasars with redshift z < 0.83 
in Data Release 3 of the Sloan Digital Sky Survey (SDSS). The sample consists of 3772 
quasars with better model of H/3 and [O III] lines and available radio luminosity, includ- 
ing 306 radio-loud quasars, 3466 radio-quiet quasars with measured radio luminosity or 
upper-limit of radio luminosity (181 radio-quiet quasars with measured radio luminosity). 
The virial supermassive black hole mass (Mbh ) is calculated from the broad H/3 line, the 
host stellar velocity dispersion (cr*) is traced by the core [O III] gaseous velocity dispersion, 
and the radio luminosity and the radio loudness are derived from the FIRST catalog. Our 
results are follows: (1) For radio-quiet quasars, we confirm that there is no obvious devia- 
tion from the Mbh — G * relation defined in inactive galaxies when Mbh uncertainties and 
luminosity bias are concerned. (2) We find that radio-loud quasars deviate much from the 
Mbh — c* relation respect to that for radio-quiet quasars. This deviation is only partly due to 
the possible cosmology evolution of the Mbh — c* relation and the luminosity bias. (3) The 

radio luminosity is proportional to M BH ~° 16 (^Boi/iEdd) 1 ' 29-0 ' 24 for radio-quiet quasars 

3 20~t~0-60 4-1 40 

and M BH "° ™ (L Bo i/L Edd ) 4 ' 18 -i -i° for radio-loud quasars. The weaker correlation of the 
radio luminosity dependence upon the mass and the Eddington ratio for radio-loud quasars 
shows that other physical effects would account for their radio luminosities, such as the black 
hole spin. 

Key words: quasars: emission lines — galaxies: nuclei — galaxies: bulges — black hole 
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1 INTRODUCTION 

The relation between the supermassive black hole (SMBH) mass and the host stellar velocity dispersion 
(hereafter Mbh — c* relation) is one of the most important results in the study of supermassive black 
holes (SMBHs) in these decades, implying the intimate correlation between the SMBHs and their host 
galaxies (e.g. Gebhardt et al. 2000; Ferrarese & Merrit 2000; Tremaine et al. 2002; Lauer et al. 2007). This 
correlation would provide strong constraints for the evolution of active galactic nuclei (AGNs) if we know 
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AGNs follow this relation or not. However it is still under debate for different kind of AGNs, such as radio- 
loud AGNs, narrow-line Seyfert 1 galaxies, intermediate supermassive black hole, et al. (e.g., Nelson 2001; 
Boroson 2003; Shield et al. 2003; Bian & Zhao 2004; Grupe & Mathur 2004; Bonning et al. 2005; Greene 
& Ho 2006; Woo et al. 2006; Zhou et al. 2006; Salviander et al. 2007; Komossa & Xu 2007; Shen et al. 
2008). In order to study this relation for AGNs, we should calculate M B h and a * as accurately as possible. 

The width of the broad emission line (e.g., H/3 , Ha , Mg II , C IV ) can be used to trace virial velocity 
of the clouds in broad line regions (BLRs) when the line contribution from narrow-line regions (NLRs) is 
reasonably removed, and the reverberation mapping method or the empirical luminosity-size relation can 
be used to calculate the BLRs size (e.g., Kaspi et al. 2000; McLure & Dunlop 2004; Bian & Zhao 2004; 
Peterson et al., 2004; Greene & Ho 2005b). The gas velocity dispersion of the narrow lines (e.g., [O III] 
> [O II] , [S II] ) from NLRs are usually used to trace the host stellar velocity dispersion (e.g., Nelson & 
Whittle 1996; Greene & Ho 2005a). We also can directly measure the host velocity dispersion from AGNs 
host spectra (e.g., Kauffmann et al. 2003; Heckman et al. 2004; Greene & Ho 2005a; Bian et al. 2006). The 
larger number of quasars found in the Sloan Digital Sky Survey (SDSS) provides the possibility to tackle 
the Mbh — cr* relation in radio-loud quasars, (e.g. Bian & Zhao 2004; Salviander et al. 2007). 

The dichotomy of the radio loudness in quasars is a long-time question since the discovery of quasars 
(Sandage 1965; Strittmatter et al. 1980; Kellermann et al. 1989). The radio luminosity is assumed com- 
ing from the relativistic electrons powered by a jet, which is intimately connected with the SMBH (e.g., 
Begelman et al. 1984; Blundell & Beasley 1998). For scale-free jet physics and accretion theories, the radio 
luminosity is related to the central engines, such the SMBH mass, the SMBH spin, the Eddington ratio, 
et al. (Heinz & Sunyaev 2003). For radio-loud or radio-quiet quasars, the dependence of the radio loud- 
ness/luminosity upon the SMBH mass/Eddington ratio is discussed by many peoples, some support it and 
some against it. (e.g., Franceschini et al. 1998; Laor 2000; Lacy et al. 2001; Ho 2002; Woo & Urry 2002; 
McLure & Jarvis 2004; Wang et al. 2004; Greene et al. 2006; Liu et al. 2006; Sikora 2007; Panessa et al. 
2007). Laor (2003) gave some comments on the origin of AGNs radio loudness and discussed the some 
error SMBH mass estimation for radio-loud AGNs in the literature, which is mainly due to optical spectra 
with low signal-to-noise ratios, no correction of H/3 contribution from narrow line regions (NLRs). 

In this paper, we use larger number of quasars with redshifts z < 0.83 in SDSS Data Release 3 (DR3; 
see Abazajian et al. 2005) to investigate the M B h — relation and the radio luminosity dependence on the 
SMBH mass and the Eddington ratio for radio-loud and radio-quiet quasars. In §2, we briefly introduce the 
SDSS quasars Data Release 3 catalog of Schneider et al. (2005). §3 is the data analysis. §4 introduces the 
methods to calculate the SMBH masses and the Eddington ratios. Our results and discussions of Mbh — 
relation and the origin of radio luminosity are given in §5 and §6, respectively. The last section is our 
conclusions. All of the cosmological calculations in this paper assume H = 70km s _1 Mpc -1 , Qm — 0.3, 
Ov = 0.7. 

2 SAMPLE AND DATA ANALYSIS 

The sample used in this paper is selected from the SDSS quasars Catalog III, which covers a spectroscopic 
area of 1360 sq. deg., about 40% of the proposed SDSS survey area (Schneider et al. 2003). This catalog 
consists of 46,420 quasars in SDSS DR3 with Mi < —22. The catalog also contains radio emission proper- 
ties from Faint Images of the Radio Sky at Twenty-cm (FIRST) survey within 2.0" of the quasars position 
(see Col. 17 in their Table 1). 

SDSS optical spectra cover the wavelength range 3800-9200 A with a resolution of 1800 < R < 2100. 
In order to calculate SMBH mass from the broad H/3 line and the host stellar velocity dispersion from 
the narrow [O III] line, we just consider the quasars with redshifts less than 0.83, which consists of 9753 
quasars. Because whether the SMBH mass from Mg II linewidth is consistent with that from H/3 line width 
is still a complex question (e.g., Salviander et al. 2007), here we don't consider using Mg II linewidth to 
calculate the SMBH mass. 

The radio luminosity at 5GHz is calculated from the peak flux density listed in Col. 17 in Table 2 
(Schneider et al. 2003), considering the spectral index of a — 0.5, where f v oc v~ a . The radio loudness 
R is calculated from: R = /sghz/Zb, where /sghz and ,/b are the rest-frame flux density at 5 GHz and 
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Fig.l Sample of SDSS spectrum measurement for SDSS Jl 13801.84+490506.5. In the top 
panel, the black curve is the observed spectrum, the red line is the sum of the power-law contin- 
uum, the Balmer continuum and Fe II multiples (blue curves). The green ranges are our fitting 
windows. The bottom panel is the multi-Gaussian fit for H/3 and [O III] lines. The red line is the 
sum of all multi-Gaussian (blue curves). The green curve is our fitting range of the pure H/3 and 
[O III] emissions after the subtraction of the power-law continuum, the Balmer continuum and 
Fe multiples. 

4400A, considering k correction. R = 10 is commonly used to define radio-loud quasars and radio-quiet 
quasars (e.g., McLure & Jarvis 2004), as well as the radio luminosity at 5 GHz (e.g., Lacy 2001). 

For 9573 quasars with z < 0.83 from SDSS DR3, 914 quasars are detected by FIRST, 7846 quasars 
are under the FIRST flux limit, and 993 quasars are not in the region covered by FIRST. For these objects 
with non-detection in FIRST, we only have the upper-limits of the radio luminosity and the radio loudness. 
598 quasars with detection in FIRST and R > 10 are classified as radio-loud quasars. 316 quasars with 
detection in FIRST and R < 10 are classified as radio-quiet quasars. 5712 quasars with non-detection in 
FIRST and R < 10 are classified as radio-quiet quasars but with upper-limits of R and the radio luminosity. 

As we know, NLRs can contribute H/3 emission in the total H/3 profile; [O III] usually shows non- 
symmetric profile and its narrow/core component can trace the stellar velocity better (e.g., Greene & Ho 
2005a); optical and ultraviolet Fe II multiples are often presented in quasars spectra; Balmer continuum is 
required because of the existence of strong Balmer emission lines, therefore, we use following steps to do 
the SDSS spectral measurements. 

(1) First, we do the Galactic extinction in the observed spectra by using the extinction law of Cardelli, 
Clayton & Mathis (1989) (IR band) and O'Donnell (1994) (optical band), then the spectra are transformed 
into the rest frame defined by the redshifts given in their FITS headers. 

(2) The optical and ultraviolet Fe II template from the prototype NLS1 I ZW 1 is used to subtract 
the Fe II emission from the spectra (Boroson & Green 1992; Vestergaard & Wilkes 2001). The I ZW 
1 template is broadened by convolving with a Gaussian of various linewidths and scaled by multiplying 
a factor. A power-law continuum and the Balmer continuum are added in the fitting. We calculate the 
Balmer continuum following Grandi (1982) and also add the high order Balmer lines at the red side of the 
Balmer edge using the result in Storey & Hummer(1995). The best subtraction of the Fe II , power-law and 
Balmer continuum is found when \ 2 minimized in the fitting windows: 3550-3645, 4170-4260, 4430-4770, 
5080-5550, 6050-6200, 6890-7010A (see a sample fit in the top panel of Figure 1). The monochromatic 
luminosity at 5100A(AL.\(5100A) ) is calculate from the power-law continuum. 

(3) Two sets of two-Gaussian are used to model [O III] AA4959, 5007 lines. Three-Gaussian is used to 
model H/3 line. For the doublet [O III] AA4959,5007, we take the same linewidth for each component, and 
fix the flux ratio of [O III] A4959 to [O III] A5007 to be 1 :3. Two components of H/3 (supposed from NLRs) 
are set to have the same linewidth of each component of [O III] A5007 and their flux are constrained to be 
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less than 1/2 of each component of [O III] A5007. The linewidth of the broad component of H/3 is used to 
trace the virial velocity around central SMBH (see a sample fit in the bottom panel of Figure 1). 

From above spectral measurement, we obtain the full width at half maximum (FWHM) of the broad 
H/3 line and the narrow/core [O III] line (FWHMh/3, FWHM" oiii j), the monochromatic luminosity at 

5100A(ALa(5100A) ), the total H/3 luminosity (Lh/jX as well as the radio luminosity and the radio loudness 
for SDSS DR3 quasars with z < 0.83. 

Objects without the H/3 or [O III] lines are eliminated. In order to obtain the reliable spectra fit, we 
carefully select objects for analysis. The line equivalent width (EW) can show line signal-to-noise ratios. 
The error of EW can be regard as a tracer to show the fitting goodness. Because the Hf3 is usually strong, 
we don't constrain EW of H/3 line, only constrain the error of EW for H/3 line. We select objects by the 
criterions of EW of [O III] larger than 1.5, the errors of EWs of H/3 and [O III] AA4959, 5007 less than 
100%. It leads to 367 radio-loud quasars, 3677 radio-quiet quasars including 207 radio-quiet quasars with 
measured radio luminosity. Then we visually check these spectra one by one. 

At last, we obtain a sample of 3772 quasars with better multi-components model of H/3 and [O III] lines, 
including 3466 radio-quiet quasars (hereafter "RQ total sample"), 306 radio-loud quasars (hereafter "RL 
sample"). Most objects in these 3466 radio-quiet quasars only have upper-limits of the radio luminosity 
and the radio loudness, 181 radio-quiet quasars (hereafter "RQ sample") have the measurements of the 
radio luminosity and the radio loudness. We use the radio-quiet sample as the control sample to discuss the 
A'Tbh — c* m radio-loud quasars. 

3 SMBH MASS, EDDINGTON RATIO AND STELLAR VELOCITY DISPERSION 

The BLRs size is calculated from the monochromatic luminosity at 5100A (\L\(5100X) ) or the H/3 
luminosity by the following formulae (Kaspi et al. 2005): 

/ \ t /rinn^W 0.69±0.05 

<f tooA) = (22 3 ± 21) lt _ days (1) 



/ j \ 0.80±0.11 

ife = (82.3 ± 7.0) It - days (2) 

We use the FWHM of the broad H/3 line (FWHM H/3 ) to trace the BLRs virial velocity u B lr = V7 x 
FWHMh/3, / is the calibration factor. If BLRs cloud is disk-like with a inclination of 8 (Wills & Browne 
1986), 

FWHMh^ = 2(« r 2 + ^LRsin^) 1 / 2 (3) 

where v r is the random isotropic component. We can then calculate the SMBH masses by Mbh = 
flBL1 ^ LR (Kaspi et al. 2000; Kaspi et al. 2005): 

/FWHM H/3 \ 2 /AL A (5100A)\ 0,6 " 



Mbh = f x 4.35 x 10° — £ f M Q . 

Vl0 3 kms-V VlO^crgs- 1 / 

(4) 



/ FWHMh/3 \ 2 f /.II-; '"' S " 



M BH = / x 1.61 x 10' ^ „ ' — T \ M Q . 



(5) 



If assuming v r <§C ?;blr, and the random orbits of BLRs clouds, / = 0.75. Onken et al. (2004) did a 
calibration by the M B h — <J* relation and suggested / w 1.4 (see also Collin et al. 2006; Dasyra et al. 
2007). In our mass calculation, we adopt the random orbits of BLRs clouds and / = 0.75. 
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Fig.2 The distributions of M B h, AL A (5100A), FWHM H/ 3, c^ni] for 181 radio-quiet quasars 
with radio loudness (top), 306 radio-loud quasars with measured radio luminosity (middle), total 
3466 radio-quiet quasars (bottom). 



We calculate the Eddington ratio, i.e., the ratio of the bolometric luminosity (Lboi) to the Eddington 
luminosity (£Edd)> where L^dd = 1-26 x 10 38 (MBH/-W©)erg s _1 . The bolometric luminosity is calculated 
from the monochromatic luminosity at 5100A , Lboi = cbAL a (5100A), where we adopt the correction 
factor c B of 9 (Kaspi et al. 2000; Marconi et al. 2004; Richards et al. 2006; Netzer & Trakhtenbrot 2007). 

We use the gas velocity dispersion of the narrow/core [O H i] component from NLRs to trace the host 
stellar velocity dispersion, Opurj ~ V a obs — [finst/(l + z)] 2 , where cr obs = FWHM™ onl ]/2.35, z is the 
redshift (Bian et al. 2006). For SDSS spectra, the mean value of instrument resolution a inst is 60 km s _1 
for [O III] (e.g. Greene & Ho 2005a). 

In Figure 2, we present the distributions of the SMBH mass, ALa(5100A), FWHM H/ a, cjbui] for 
181 radio-quiet quasars with radio loudness (top), 306 radio-loud quasars with measured radio luminosity 
(middle), total 3466 radio-quiet quasars (bottom). The mean of SMBH mass is 8.65 ± 0.03 with a standard 
deviation of 0.45 for RL sample of 306 radio-loud quasars, 8.36 ±0.04 with a standard deviation of 0.48 for 
RQ sample of 181 radio-quiet quasars with reliable radio luminosity, 8.32 ± 0.01 with a standard deviation 
of 0.43 for total 3466 radio-quiet quasars. Radio-loud quasars have larger SMBH masses, and there is only 
a few objects with mass less than 10 8 M© (see Figure 3), which is consistent with the results of McLure & 
Jarvis (2004). Radio-loud quasars have smaller Eddington ratios, respect to radio-quiet quasars (see Table 
1). We find that, for radio-loud quasars, the mean of H/3 FWHM is 7493 ± 165 km s _1 with a standard 
deviation of 2882 km s^ 1 , the mean of log AL A (5100A) is 44.86 ± 0.03 erg s" 1 with a standard deviation 
of 0.45; for radio-quiet quasars, the mean of H/3 FWHM is 5780 ± 176 km s _1 with a standard deviation 
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Fig. 3 The Mbh — c* relation for radio-loud and radio-quiet quasars. Red circle denotes radio- 
loud quasars, blue star denotes radio-quiet quasars with measured radio luminosity, black square 
denotes the radio-quiet quasars with upper-limit of the radio luminosity. The mass in the left panel 
is derived from AL>(5100A) , and the mass in the right panel is derived form Hf3 luminosity. 

of 2389 km s _1 , the mean of log ALa(5100A) is 44.81 ± 0.01 erg s _1 with a standard deviation of 
0.46. Radio-loud quasars tend to have larger H/3 FWHM and ALa(5100A) , leading larger SMBH masses 
(Sulentic et al. 2000). 

4 Mbh ~ RELATION 

4.1 The mass deviation from Af B H — cr* relation 

In Figure 3, we show Mbh — &* relation for radio-loud and radio-quiet quasars. The solid line in Figure 
3 is the A/bh — a* relation in normal nearby galaxies given by Tremaine et al. (2002), Mbh(c*) = 
10 8 - 13 [cr*/(200 kms" 1 )] 4 02 M Q . In the left and right panels of Figure 3, AL A (5100A) and L Hf3 are used 
to calculate the SMBH mass, respectively. In Figure 3, the correlation between A/bh and <J™ OUVi is very 
weak for larger SDSS quasars sample. It is possibly due to the accuracy of the stellar velocity dispersion 
derived from the narrow/core [O III] line- width. However, it is obvious that the sample of radio-loud quasars 
deviated much from the solid line respect to that for the sample of radio-quit quasars. It is consistent with 
our previous result (Bian & Zhao 2004). 

We calculate the black hole mass deviation AlogA/BH from the solid line defined by Tremaine et al. 
(2002), AlogA/BH = logAiBH(H/3) — logMBH(c*), where cr* is adopted to be ajoiiil- For the mass derived 
from ALa(5100A) , the mean of AlogA/BH 1 is 0.65 ±0.04 with a standard deviation of 0.71 for RL sample 
of 306 radio-loud quasars, 0.04 ± 0.04 with a standard deviation of 0.63 for RQ sample of 181 radio-quiet 
quasars with reliable radio luminosity, 0.14 ± 0.01 with a standard deviation of 0.62 for total 3466 radio- 
quiet quasars. We find that they are almost the same for the case of the mass derived from L-^p- In the next 
analysis, we just consider the mass and Eddington ratio calculated from XL\ (5100A) . 

In the top panel of Figure 4, we plot the deviation of the SMBH mass from Mbh — c* relation versus the 
radio loudness. It is obvious that the deviation tends to be larger when radio loudness becomes larger. In the 
bottom panel of Figure 4, we also plot the deviation of the SMBH mass from the Mbh — c* relation versus 
the redshift. We find a weak correlation between the mass deviation and the redshift for radio-quiet quasars. 
The simple least-square regression gives: AlogA/BH = (1-00 ± 0.06)z — (0.29 ± 0.03). The correlation 
coefficient R is 0.26, with a probability of p n ull < 10~ 4 for rejecting the null hypothesis of no correlation. In 
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Fig. 4 Top: The deviation of the SMBH mass from the Tremaine's Mbh — o* relation in fig 2. 
versus the radio loudness. The dash line denotes AlogAfsH = 0. Bottom: The deviation of the 
SMBH mass from the Tremaine's Mbh — o - * relation in fig 2. versus the redshift. The red solid 
line denotes our best fit for all radio-quiet quasars. Symbols as Figure 3. 



Figure 5, we show the redshift distributions for radio-quiet and radio-loud quasars. The radio-loud quasars 
(red circles) have larger redshifts relative to the radio-quiet quasars (blue stars) (see Figure 5). 

In Table 1, we show the mean values of the masses and the Eddington ratios in different redshift bins 
for different samples. 

4.2 Uncertainties 

There are some factors to account for the uncertainty of the SMBH mass calculation: the uncertainties 
of H/3 , [O III] line widthes, ALa(5100A) , L^p when the multi-components are used to model SDSS 
spectra; the system errors in equations (1-5) from the uncertainties of the BLRs geometry and dynamics. 
The uncertainty of our calculated SMBH mass is about 0.5 dex. The uncertainty of the Eddington ratio is 
about 0.5 dex or more. For radio-loud quasars, we should account two effect: the relativistic beaming effect 
on the optical continuum and the orientation of BLRs. The total H/3 luminosity instead of \L\(5100X) 
is used to account for the first effect. We find the effect is small in our sample, and there is no correlation 
between the H/3 EW and the radio loudness (e.g. Wu et al. 2004). Lacy et al. (2001) made a small correction 
of the orientation of BLRs by a factor of R® 1 , where R c is the ratio of core to extended radio luminosity. 
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Fig. 5 The redshift distributions for 1 8 1 radio-quiet quasars with radio loudness (top), 306 radio- 
loud quasars with measured radio luminosity (middle), total 3466 radio-quiet quasars (bottom). 



They adopted R c = 0.1 for steep- spectrum quasars and R c = 10 for flat-spectrum quasars if R c is not 
measured. It will lead to the uncertainties of AlogMBH about 0.2 dex. 

The fibers in the SDSS spectroscopic survey have a diameter of 3" on the sky. The SDSS spectra 
of lower-redshift quasars possibly have obvious stellar light contribution, which can be used to directly 
measure the stellar velocity dispersion (e.g., Kauffmann et al. 2003; Heckman et al. 2004; Bian et al. 2006). 
For luminous quasars (Mi < —22), the stellar light contribution can be omitted or has little effect on the 
mass calculation (e.g., Vanden Berk et al. 2006). 

It is possible that jet can have a dynamical effect on the NLRs and may have a systematically different 
effect on the [O III] profile (Nelson & Whittle 1996). However, considering the [O III] profile broadening 
by jet, the correction of [O III] gas velocity dispersion will lead the radio-loud quasars to deviate much 
more from the Mbh — er* relation. 

We select quasars with the EW of narrow [O III] component larger than 1.5, EW errors of H/3 and [O 
III] AA4959, 5007 less than 100%. Different criteria would lead to different number of quasars (such as 
error of EWs less than 5%, 50%, or 100%, x 2 < 4). However, we find that the main results don't change. 
In the next subsection, we will discuss the luminosity bias in detail. 

4.3 The mass deviation from the luminosity bias 

Salviander et al. (2007) also used SDSS DR3 quasars to discuss the cosmological evolution of Mbh — 
er« relation. After carefully consider the selection biases and intrinsic scatter in the Mbh — cr* relation, 
they suggested that Mbh — cr* relation appears to evolve with redshift. Netzer & Trakhtenbrot (2007) 
also found the nonlinear Mbh — relation with the different slopes for different redshift bins. During 
our process of selecting objects, the line fitting favor the brighter objects (i.e. luminosity bias, see Figure 
6). Following the work of Salviander et al. (2007), we calculate the contribution of AlogMBH from this 
luminosity bias. We calculate the mean observed luminosity in different redshift bins for our different 
samples (i.e. RQ total sample; RQ sample; RL sample). Using the QSOs luminosity function (Boyle et 
al., 2000), we calculate the cut luminosity to make the mean luminosity of the kept QSOs (L > L cut ) 
equal to the observed mean luminosity in different redshift bin. Then we do the simulation to calculate the 
contribution of AlogMBH from this luminosity bias (for detail in Salviander et al. 2007). We obtained a 
formulae: AlogMff u = 0.292 + 0.1138a; + 0.265a; 2 + 0.480a; 3 + 0.182a; 4 , where x is log{L cut /L ), 
L Q = 0.3L E dd(M* ), M* j = lO n M (Drory et al. 2005). We find that the mass deviation form the 
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Fig. 6 ALa(5100A) versus z. Open circles denote RL sample; blue stars denote RQ QSOs with 
measured radio loudness; blue squares denote RQ QSOs with upper-limits of radio loudness; 
small black squares denote all 9753 SDSS DR3 QSOs with z < 0.83. It is obvious that some 
faint objects are missed in our selection. 
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0.1-0.4 


114 


44.35 ± 0.32 


8.21 ±0.32 


-1.01 ±0.39 


-0.04 ±0.63 


-1.14 


0.10 


0.4-0.82 
RL 


67 


44.94 ± 0.34 


8.62 ±0.41 


-0.83 ±0.32 


0.17 ±0.61 


-0.50 


0.25 


0.1-0.3 


35 


44.25 ± 0.29 


8.25 ±0.38 


-1.15 ±0.42 


0.01 ±0.68 


-1.22 


0.08 


0.3-0.5 


109 


44.49 ± 0.31 


8.51 ±0.45 


-1.17 ±0.39 


0.52 ± 0.80 


-0.98 


0.15 


0.5-0.7 


106 


44.78 ± 0.35 


8.78 ±0.38 


-1.14 ±0.32 


0.80 ± 0.63 


-0.66 


0.23 


0.7-0.83 


56 


44.98 ± 0.29 


8.91 ±0.39 


-1.08 ±0.33 


0.95 ± 0.75 


-0.46 


0.26 



Table 1 The mean quantities in different redshift bins for different samples. L cut is calcu- 
lated from the QSOs luminosity function (Boyle et al. 2000) to make the mean luminosity 
of the kept QSOs (L > L cut ) equal to the observed mean luminosity in different redshift 
bins. L Q = 0.3LEdd(-Mg a i), where M* al = 10 11 M Q in the galaxy mass function $(A/ ga i) = 
$*(M gal /Af* al )" a e~ M ^'/ M g *ai (Drory etal. 2005). 



luminosity bias is monotonously increased with the redshift (see Col. (8) in Table 2). Table 1 shows our 
results. Col.(l) is the redhsift bin; Col. (2) is the number in the redshift bin; Col.(3)-(5) are the mean values 
of 5100Aluminosity, mass, Eddington ratio; Col. (6) is the mean mass deviation from Mbh — cr* relation; 
Col. (7) is log(L cut /Lo); Col. (8) is our simulated mass deviation for different cut luminosity in different 
redshift bins. 
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Fig. 7 The radio luminosity versus the [O III] luminosity. Red circle denotes radio-loud quasars, 
blue star denotes radio-quiet quasars with measured radio luminosity, black square denotes the 
radio-quiet quasars with upper-limit of the radio luminosity. The red solid line is the BCES 
bisector result for radio-quiet quasars (blue stars). The blue solid line is the BCES bisector result 
for radio-loud quasars (red circles). The dash line is the best fit of radio-quiet AGNs found by 
Xuetal. (1999). 



For RQ total sample, the luminosity bias can interpret most amount of AlogAfeH- For the highest 
redshift bin of 0.7 < z < 0.83, AlogA/sH is about 0.15 dex after correction the effect of luminosity bias. 
This 0.15 dex is possibly the A/bh — o* cosmological evolution in this highest redshift bin, which is very 
consistent with the result of Salviander et al. (2007). However, we should note that the standard deviation 
of AlogA/eH in different bins is about 0.6dex, which is very larger than this 0.15dex. For RQ sample, the 
observed AlogAfeH can be completed contributed from the luminosity bias, which is possibly due to the 
smaller numbers of this sample. Therefore, we think there is no obvious deviation from A/bh — relation 
considering the A/bh uncertainties and the luminosity bias. 

For RL sample, after corrected the contribution from the luminosity bias, A log A/bh is still large 
(about 0.69 dex in 0.7 < z < 0.83) and there is a trend that AlogAfsH becomes larger for larger redshift 
bin. Considering the possible Mbh — u* cosmological evolution (0.15 dex in 0.7 < z < 0.83), for radio 
loud QSOs, there are still 0.54 dex deviation in 0.7 < z < 0.83. Bonning et al. (2005) suggested that 
narrower [O III] for radio loud quasars is responsible for this deviation from the A/bh — 0* relationship, 
and it is not the effect involving Mbh ■ The cause of this deviation is unclear. 

5 ORIGIN OF RADIO LUMINOSITY 

5.1 L 5G Hz - £[oni] relation 

The relation between the radio luminosity and the optical/X-ray luminosity, which provides the connection 
between the jet and accretion power, have been discussed by many group (e.g. Xu et al. 1999; Ho 2002; 
Wang et al. 2004; Panessa et al. 2007; Sikora 2007). In Figure 7, we show the radio luminosity at 5GHz 
versus the total [O III] luminosity. These two luminosities are all related to the redshift. By the partial 
Kendall's r correlation test, we do a partial correlation analysis with redshift as the test variable(Akritas 
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& Siebert 1996). For the RQ sample of 181 radio-quiet quasars, partial Kendall's r correlation is 0.237, 
variance is 0.0304, and the probability of null hypothesis is 6.3 x 10~ 15 . For RL sample of 306 radio-loud 
quasars, r correlation is 0.251, variance is 0.0456, and the probability of null hypothesis is 3.7 x 10~ 8 . We 
use the bivariate correlated errors and intrinsic scatter (BCES) regression methodQof Akritas & Bershady 
(1996) (see also Isobe et al. 1990) to find the relation between L\q ln j and £501^ and adopt the BCES 
bisector result (e.g. Kaspi et al. 2005). For RQ sample of 181 radio-quiet quasars with measured radio 
luminosity, the BCES bisector result: logi 5G Hz = (0.87±0.04)logL [o H i] + (3.40± 1.87) (red dash line in 
Figure 7). ForRL sample of 306 radio-loud quasars, logL 5 GHz = (1.60±0.08)logL[ O m] — (26.30 ±3.39). 

In Figure 7, considering the errors of the intercept, our best fits for radio-quiet quasars is consistent with 
the result found by Xu et al. (1999) (also see Ho & Peng 2001) : logL 5G Hz = (0.45 ± 0.07)logL [o in] + 
(20.25 ± 0.6) (black dot line in Figure 7). In the plot of radio luminosity versus the optical/X-ray nuclear 
luminosity, the separation of radio-loud and radio-quiet quasars from SDSS DR3 is not too clear as other's 
results (Xu et al. 1999; Terashima & Wilson 2003; Sikora et al. 2007). The difference is possibly due to the 
selection effect by different wavelength bands. 

[O III] luminosity is usually assumed to be proportional to the accretion rate and this correlation can be 
explained in a model of accelerated and collimated jet by magnetic field (Xu et al. 1999). Apart from the 
dependence on the accretion rate, the radio luminosity possibly dependents on the central SMBH properties, 
mass or spin (e.g., Sokira et al. 2007 and refs. therein), which we will discuss in the next section. 

If we use the tight correlation between X-ray luminosity and [O in] luminosity (e.g., Xu et al. 1999), 
logL x = 1.011og\L[om] + 1.6, the relation between Lsghz and L[om] can be transformed to the relation 
between L 5G h z and L x : L 5G Hz oc £0.86±o.06 f or RQ samp ie and £ 5G Hz oc L x - 58±01 ° for RL sample. There 
exists obvious different on the slope for radio-quiet and radio-loud quasars. For low luminosity AGNs, 
Panessa et al. (2007) suggested a correlation, L x oc L®qJ^ z , their index is between ours for radio-quiet and 
radio-loud quasars. If we use the correlation suggested by Netzer et al. (2006), £oin oc £0- 704 ±°- 06 5 the 
relation we found between £sghz and L[om] can be transformed to -Lsghz oc £°- 61±0 04 for RQ sample 
and £5 G Hz oc ^ 11±0 07 for RL sample. The radio luminosity is often assumed coming from the relativistic 
electrons powered by a jet. The result of RL sample is consistent with that of Panessa et al. (2007). The 
X-ray emission is often assumed coming from both the accretion flow and the relativistic jet, dominated by 
accretion flow at high accretion rate, and dominated by jet emission at low accretion rate. (Gallo et al. 2003; 
Yuan & Cui 2005). This relation between Lsghz and L x in different accretion rates can be explained in the 
jet-dominant X-ray models (e.g., Fender et al. 2003; Gallo et al. 2003; Heinz 2004; Yuan & Cui 2005). 

5.2 The radio luminosity dependence on the SMBH mass and the Eddington ratio 

It is suggested that the radio luminosity/radio loudness is relate to the SMBH masses (e.g., Laor 2000). We 
calculate the radio luminosity dependence on the SMBH mass and the Eddington ratio, i.e. logLsGHz = 
a + 6ilogM B H + ^2log(L Bo i/iEdd) (see Figure 8). 

We firstly do the multiple regression with ASURV Rev 1.2 (LaValley, Isobe & Feigelson 1992 and refs. 
therein) for RQ sample, RL sample and RL+RQ sample. In order to avoid the non-symmetric regression, 
a,bi,b2 are adopted the mean values when we use different variable as the dependent variables in the 
multiple regressions (see Col.(l) in Table 1). In all the multiple regressions, the probability for rejecting the 
null hypothesis of no correlation is p nu \\ < 10~ 4 . The R-Square correlation coefficient for RQ sample is 
larger than other two samples (see Table 2, Figure 8). 

We also do the symmetric multivariate regression analysis, through the equation y = a + bixi + 62^2, 

directly by the X 2 estimator, X 2 = E 4 J ^'^y^b^ )2 ' (Press et aL 1992; Tremaine et al - 2002 ; 
Merloni et al., 2003), where a are the corresponding uncertainties. Considering the same uncertainties (a) 
of radio luminosity, mass, and the Eddington ratio (Tremaine et al., 2002; Merloni et al., 2003), we re- 
normalized these uncertainties to make the minimum x 2 / n doi of unity, the results are listed in first lines for 
different sample in Table 1 . 

Considering the errors of a, b%, 62 in Tables 2, the results from ASURV and \ 2 are consistent very well. 
Therefore, in the next analysis, we adopt the values of a, bi, 62 from x 2 estimator (Table 1), i.e. Lsghz oc 



This is not a symmetric regression used by Merloni et al. 2003. For detail in section 5.2 
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Fig. 8 The radio luminosity dependence on the SMBH mass and the Eddington ratio. The indexes 
are adopted from the mean values in brackets in Table 1 . Top panel is for radio-loud quasars, 
middle panel is for radio-quiet quasars with measured radio luminosity, bottom panel is for them 
all. Symbols as Figure 7. The dash lines denote 1:1. 



1 28 +0 31 3 10 +1 40 

M BH - Q ^OWiEdd) 1 - 29 - ^ for RQ sample, and L 5GHz <x Af BH "° 70 OW^Edd) 4 ' 18 - 1 10 for RL 

4 30 +0 ' 70 c i ,- + 2.32 

sample, and L 5G h z oc M bh "° 80 (L B oi/L E dd) 1 m for RL+RQ sample. 

Ho (2002) suggested a correlation between the nuclei radio loudness and the Eddington ratio (Gallo et 
al. 2003; Greene et al. 2006; Sikora et al. 2007; Panessa et al. 2007). We also use the multiple regression 
by ASURV to search the radio loudness dependence on the SMBH mass and the Eddington ratio (set radio 
luminosity as the dependent variable). However, the R-Square correlation coefficient is very low for RL 
sample and RQ sample. For the RL+RQ sample, we find a weak correlation between the radio loudness 
and the SMBH mass (the simple least-square correlation R=0.26), much weaker correlation between the 
radio loudness and the Eddington ratio (R=-0.15). The range of Eddington ratio is between 0.01 to 1 for 
our RL sample and RQ sample. And our sample is composed by broad line type I quasars, which just fill 
the gaps between two sequences in plot of radio loudness versus the Eddington ratio (see Figure 3 in Sikora 
et al. 2007). When we research the disk-jet connection model, X-ray luminosity is maybe a better tracer of 
SMBHs accretion power than the optical luminosity (e.g. Panessa et al. 2007). We also should pay more 
attention on narrow-line Seyfert 1 galaxies with larger Eddington ratios in this kind of plot (e.g. Zhou & 
Wang 2002; Whalen et al. 2006; Komossa et al. 2006). 
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Table 2 The mean quantities in different redshift bins for different samples, a, b±, 62 is defined 
by: logL 5 GHz = a + 6ilogMeH + folog^Boi/^Edd)- For different samples, the first line is for 
the result by \ 2 minimization, and in the second line, quantities in brackets are the mean values 
of a, bi, 62 by ASURV, the last three lines are results considering different dependent variable by 
ASURV. 



Dependent variable 


a 


bi 


bi 


R-sqaure 


(i) 


(2) 


(3) 


(4) 


(5) 


RQ 


30.9±i;* 


J-^°_0.16 


1 29 +u ' 31 






(30.84 ±2.21) 


(1.27 ±0.29) 


(1.26 ±0.53) 




I/5GHz 


33.18 ±0.43 


0.95 ± 0.05 


0.81 ±0.07 


0.64 


M BH 


28.79 ± 1.54 


1.52 ± 0.02 


1.12 ±0.08 


0.72 


Z/Bol /-t/Edd 


30.54 ± 3.04 


1.35 ± 0.09 


1.85 ±0.02 


0.57 


RL 


ly - ' -3.90 


3 10+ ,) ■ ^i,l 

°- ±u -0.70 


4.18i{;3 






(22.43 ± 10.4) 


(2.68 ± 1.30) 


(3.29 ± 1.94) 




£5GHz 


32.62 ±0.82 


1.24 ±0.10 


1.38 ±0.13 


0.35 




11.77 ±3.54 


3.85 ± 0.40 


3.23 ±0.19 


0.53 


Z/Bol /-t/Edd 


22.89 ± 3.79 


2.95 ±0.16 


5.26 ± 0.05 


0.49 


RL+RQ 


10 o +8 - YU 

1U.U_ 4 20 


4 30 +,);/u 


c 1 z+'2.-i'2 
d - id -1.69 






(15.65 ± 12.93) 


(3.46 ± 1.81) 


(4.11 ±3.81) 




L$GHz 


30.45 ± 0.73 


1.38 ±0.09 


0.92 ±0.12 


0.31 


M BH 


6.57 ±2.74 


4.34 ± 0.04 


3.09 ±0.17 


0.53 


I/Bol /LEdd 


9.92 ±4.67 


4.67 ±0.25 


8.33 ±0.08 


0.40 



For scale-free jet physics, Heinz & Sunyaev (2003) derived the dependence of the accretion-powered 
jets flux (f v ) upon the SMBH mass and the dimensionless accretion rates for different accretion sce- 
narios (see their Table 1). For radiation-pressure-supported standard accretion disk, /„ oc Mgg 12 ""'' 3 ; 
for gas-pressure-supported standard accretion disk, f v oc ^^ 7 ~ 32Q )/ 120 ? 7 l ( 1 '7i2+2a/3)4/5. f or ADAF, 
f v oc Mg^ 12 "/ 3 7 |j 17 / 12 + 2q / 3 5 where a is the radio spectral index. Assuming a = 0.5, for radiation- 
pressure-supported standard accretion disk, f v oc Mg^ 5 ; for gas-pressure-supported standard accretion 
disk, f v oc Mgjf 3 ™ 1,40 ; for ADAF, f v oc Mj^ 5 ™ 1,75 . Considering large scatter in bi,b 2 , our results are 
consistent with above the radio origin of scale-free jet model. However, by our data, we can't distinguish 
the different disks for radio-quiet and radio-loud quasars in this accretion-powered jet model. 

6 CONCLUSIONS 

With the large number of quasars in SDSS DR3 catalog, we use the multi-components to model the SDSS 
spectra and calculate the SMBH masses. Combined with the radio properties from FIRST, we obtained a 
sample of 3772 quasars with reliable SMBH masses, including 306 radio-loud quasars, 3466 radio-quiet 
quasars with measured radio luminosity or upper-limit of radio luminosity (181 radio-quiet quasars with 
measured radio luminosity). Two main results are suggested: (1) The radio-loud quasars deviate much from 
the Mbh — <?* relation of the nearby normal galaxies defined by Tremaine et al. (2002) respect to that for 
radio-quiet quasars, which is only partly due to the possible cosmology evolution of the Mbh — cr* relation 
and the luminosity bias. (2) The radio luminosity is correlated to the central SMBH mass and the Eddington 

ratio, oc Af BH ~° 16 (LBoi/^Edd) 1 ' 29 - 24 for radio-quiet quasars and oc M BH "° ™ (LBoi/^Edd) 4 ' 18 - 1 10 
for radio-loud quasars. Weaker correlation coefficient of the radio luminosity dependence upon the mass 
and the Eddington ratio for radio-loud quasars shows other physical effects would account for their radio 
luminosity, such as the SMBH spin. 
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